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The generalized Gibbs ensemble has been shown to be relevant in the relaxation of a completely 
integrable system subject to a quantum quench, in the sense that it accurately predicts the steady 
values of some physical variables. We proceed to further question its relevance by giving the quenched 
system a second quench. The concern is whether the generalized Gibbs ensemble can also accurately 
predict the relaxed system's response to the second quench. Two case studies with the transverse 
Ising model and the hard-core bosons in one dimension yield an affirmative answer. The relevance 
of the generalized Gibbs ensemble in the non-equilibrium dynamics of integrable systems is then 
greatly strengthened. 
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Recently, non-equilibrium dynamics of many-body sys- 
tems has attracted a lot of attention One common 
concern is whether an initially out-of-equilibrium system 
can thcrmalize to behave like a textbook Gibbs ensem- 
ble, and how integrability 0, Q or non-integrability of 
the system will affect its relaxation dynamics. An im- 
portant achievement on this issue is identification of the 
relevance of the generalized Gibbs ensemble (GGE) in 
the relaxation dynamics of a completely integrable sys- 
tem 01 . The so called generalized Gibbs ensemble is con- 
structed according to the principle of maximum entropy 
Q while taking into account all the constants of motion, 
whose values are determined by the initial state. With 
the same philosophy behind the construction, it is a nat- 
ural counterpart of the usual Gibbs ensembles for a non- 
intcgrable system. So far, the GGE has been found to 
predict correctly the asymptotic values of physical vari- 
ables in a variety of integrable systems [l|, |3, ■ 

The fact that asymptotically, the true, constantly 
evolving system agrees well with the GGE on the phys- 
ical quantities is definitely a non-trivial and pleasant 
one. However, one should not be content with this fact 
only. Our daily experience in the (mostly non-integrable) 
macroscopic world is that, if a system relaxes to some 
steady state, it relaxes in the sense that not only its 
static properties (i.e. values of the physical quantities) 
but also its dynamical properties agree with the steady 
state. To be specific, the system should respond to later 
perturbations as if it were indeed in the steady state. 
Therefore, it is necessary to check whether the GGE has 
this merit. If so, it surely adds to the relevance of the 
GGE in the non-equilibrium dynamics of a completely 
integrable system. It would mean that the true system 
is hardly distinguishable from the GGE neither by static 
nor dynamical criterions, and it would be fair to say the 
system has thermalized as much as possible. 

Motivated by this problem, we have studied the trans- 
verse Ising model and the hard-core bosons in one- 
dimension (which can also be mapped to the XX model) 



individually. The two models are integrable and both 
have been shown to admit a GGE account of their asymp- 
totic behaviors after a quantum quench. Here our idea 
is to give them a second quench when they have reached 
the steady phase The concern is whether they will 
respond as if the systems were in the GGE states. The 
result turns out to be the case. 

Transverse Ising model. — The Hamiltonian of the 
model is H{g) = — J2h=i i'^f'^f+i ~ d'^i)^ where af'^ are 
Pauli matrices acting on a 1/2-spin at site I. Here pe- 
riodic boundary condition is assumed and N is an even 
integer large enough. Below, quenches of the system cor- 
respond to changing the value of g (strength of the trans- 
verse magnetic field) suddenly. We will consider a double 
quench scenario. Initially the value of g is 170 and the sys- 
tem is in its ground state |Go). Then the value of g is 
changed successively to gi and 172 • 

Under the Jordan- Wigner transform {af + 

i(jf,a-f - icrf)/2 = {al,ai)exp{-iTrJ2lz}^alar), 
af = 2ajai — 1, where aj and a; are fermionic 
operators, the Hamiltonian is rewritten as 
H{g) = -J2^^^[{a+a'l^^+alai+i+h.c.)-2galai], 
with a constant term dropped [1^. Note that here 
the boundary condition is anti-periodic [l6| . Taking 
the Fourier transform bk = -^^i^^'^^^''^'^ o^U with 
k = -Ar/2-f 1/2,--- ,-l/2,l/2,-- - ,N/2 - 1/2 so as to 
comply with the anti-periodic boundary condition, we 
can rewrite the Hamiltonian as {ipk = iirk/N) 

^id) = X! [^(-9 " cos(j)k)blbk - is\i\(t>k{b^khk + ^Lfc^*!) 

k 

It is ready to verify that bk and ^L^. are coupled in 
their equations of motion and this suggests the Bo- 
goliubov transformation rjk = Ukbk + ivkb''_k- With 
Efc = 2a/H-52 _ 2gcos(j)k > 0, (cos6'fc, sin6'fc) = 2{g - 
cos0fc,sin(/)fc)/efc, and {uk,Vk) = (cos%,sin^), the 
Hamiltonian is finally diagonalized as H = J2k ^kvlvk- 
Here again the constant term is dropped. Note that Uk, 



2 



Vk, 9k, and Sk all depend on g. The dependence will be 
displayed explicitly when necessary. 

We are interested in the correlation functions (o-fcrj), 
(afaj), and the transverse magnetization (Mz) = 

ij^i'^i) — (J^kC^^k^k — !))• Here the expectation val- 
ues may be taken with respect to various states as shown 



below. Introducing Ai 



- ai and Bi 



ai, we can 



rewrite them as (afa-j) — {BiAi+iBi+i ■ ■ ■ Aj-iBj-iAj) 
and {cTiCrj) — (BiAiBjAj) [15;]. These forms allow us to 
use Wick's theorem to do the calculation. The correla- 
tion functions will be decomposed into sums of products 
of the basic correlators (AiAm), {BiBm), and {BiA,n). 

The initial state \Go) is defined as rik{ga)\Go) — for 
all k, or explicitly, [Go) oc Ofc '?fc(ffo)|V') where \tp) can 
be an arbitrary state as long as T?fc(5o)|V') 0- After 
the first quench of changing g from go to gi at t = 0, 
we have {Go\Ai{t)Am{t)\Go) Sim for t large enough 
Hill, and similarly {GQ\Bi{t)Bm{t)\Go) Sim for 



t large enough. But {Go\Bi{t)Amit)\Go) Gj 
has the value of 



(1) 



- TV ^ 



1 



k{gi) 



cos(A6i^° 



(1) 



then indicates that the asymptotic values of the corre- 
lation functions ^ and ([3]) can be recovered with the 
GGE. Likewise, the asymptotic value of the transverse 
magnetization (|4]) is exactly predicted by the GGE, i.e., 
Mi^^ =tr{MzPggei) UM- 

Now consider giving the system a second quench, 
i.e., changing the value of g from gi to 172 at some 
time t = ti. It is tedious but straightforward to 
show that at the time of t — ti + t2, for large t2 
[l6| . (Go\Ai{t)Amit)\Go) ~ (5;,„+osciUating terms de- 
pending on ti, and similarly {Go\Bi{t)Bm{t)\Go) ~ 
— <5im+oscillating terms depending on ti. However, 
{Go\Bi{t)Am{t)\Go) ~ G; ,^-l-oscillating terms depend- 
ing on ii, where 

k 



which As for the transverse magnetization, {Go\Mz{t)\Go) 



(2) 

Mz -|-oscillating terms depending on ti, with 



k 



cos0fc(52)cos(A0f )cos(A0f). (7) 



Here and hereafter 



6k{gi) - Okigj)- Thus for t 
large enough, (Gokf (t)aJ(t)|Go) ^ Gf^-: 



Gfj = Dot 



(1) 



G 



(1) 



iA+2 
^(1) 



G 



(1) 



1 ^(1) ^yi-) 



(1) 



(2) 
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and 



(1)^(1) 



{Go\a!^)a](t)\Go)^C^^=G\^>G)] 



(3) 



As for the transverse magnetization, {'i^Q\Mz{t)\'i>Q) has 
the asymptotic value of 



k 



cos0fc(gi)cos(A0f) 



(4) 



On the other hand, from go to 171, the (first) GGE density 
matrix is defined as 



Pag 



exp 



->^^k^vU9i)Vki9i)) , (5) 



with the Lagrange multiplier A^,^-* determined by the 
condition {Go\ill{gi)Vk{9i)\Go) = tr{r]l{gi)r]k{gi)pggei), 



and Zi 



n,(i 



It can be verified that 



{AlAm)ggel — —{BlBm)ggel — ^im, and {BiAm)ggel — 

g|^2i- Here the subscript means averaging over Pggei- 
Thus the basic correlators are of the same values with 
respect to the GGE density matrix Pggei and the evolv- 
ing state e~*^'^i''*|Go) for t large enough. This fact 



The oscillating terms depending on ti consist of 0{N) 
components of different non-zero frequencies and thus 
they virtually vanish for ti large enough. There- 
fore, for ti and t2 large enough, the correlation func- 
tions {Go\a'fit)a]it)\Go) and (Go|crf (t)cr|(t)|Go) have 

the same form as Eqs. ^ and ([3]) but with G^\ replaced 

by G^^]i, and {Go\Mz{t)\GQ) has the value of M^^l 

On the other hand, if the second quench is im- 
posed on the first GGE density matrix Pggei, we have 
the same asymptotic behaviors of the basic correla- 
tors and the transverse magnetization for large ^2- 

That is, {Aiit2)Amit2)) = -{Bi{t2)Bm{t2)) ^ Sim, 

{Bm{t2)Ai{t2)) ^ G^^„ and {Mz{t2)) ^ mP 0. Here 
iAiit2),Biit2),Mzit2)) = e»^(9^)*^(Ai,i?i,M,)e-*^(9^)*= 
and the average is taken over Pggei- We see that the 
transverse magnetization as well as the basic correlators 
possess the same asymptotic values regardless of the ini- 
tial state being e~*^^^i)*i |Go) or pggei- The latter fact 
implies that the correlation functions have the same prop- 
erty. However, it is not only the asymptotic values that 
can be accurately reproduced by using Pggei as a substi- 
tute for e~*^^^i-'*i |Go}. In Fig.[Tl the transient dynamics 
of Mz after the second quench is shown. There we see 
that as long as ti is large enough, the relaxation dynam- 
ics of Mz (the correlation functions have the same prop- 
erty; see the supplementary material) is independent of 
ti and can be reproduced by Pggei even to minute details. 
Therefore, as long as ti is large enough, or as long as the 
second quench comes when the system has equilibrated 
to agree with the first GGE Pggei after the first quench, 
the model reacts as if it were indeed in the GGE state 



3 




FIG. 1: (Color online) Evolution of the transverse magne- 
tization Mz after the second quench. The parameters are 
(iV, go, ffi, 32) = (10000, 2, 1, 0.2). All the lines, with the "ini- 
tial" state being the (first) generalized Gibbs ensemble (GGE) 
density matrix Pggei or e"'^'"^'*^ IGo), collapse into one. 
Here the values of t\ are chosen randomly from [500, 2500] . 
The horizontal dotted line indicates the predicted asymptotic 
value ([7|. The insert shows the time evolution of after 
the first quench. 

Pggei- That is, the GGE density matrix Pggei can serve 
as a pseudo-initial state to the second quench. 

Finally, for the quench of Pggei, we can define a second 
GGE density matrix as 

Pgge2 = yW^^^ [->^'k^vl{92)Vk{92)j , (8) 

^ A; 

(2) 

with the parameter Aj, determined by the condition 
tr{r]l{g2)'nk{g2)Pgge2) = tr{ril{g2)r]k{g2) Pggei) , and Z2 = 
+ e~ k ). The point is that the basic correlator 

Gi ^ in ([5]) and the transverse magnetization in ^ can 
be exactly reproduced by Pgge2 ■ This is one more support 
of the argument that Pggei can serve as a pseudo-initial 
state to the second quench. 

Expansion of hard-core bosons in a one dimensional 
lattice. — To make contact with previous works, the sce- 
nario studied below is an extension of that in Ref. [1]. 
There are N hard-core bosons and there is a lattice of 
M2 sites, which are numbered from 1 to M2. Initially the 
N bosons are confined to the Mq middle sites by hard- 
walls on the two sides and the system is in the ground 
state, which is denoted a.s i/jq. At t = 0, the hard- walls 
are suddenly moved outward symmetrically so that now 
Ml sites are contained. The system then evolves and as 
found by Rigol et al. the GGE plays an important 
role in the ensuing dynamics — the momentum distribu- 
tion of the bosons in its steady value is accurately cap- 
tured by the GGE density matrix S^gei (see below). Our 
idea is then at some time ti, when the momentum distri- 
bution has settled down to its steady value, to increase 
the volume to M2 sites and let the bosons expand once 
again. The aim is to see whether the subsequent dynam- 
ics can be accurately reproduced with the initial state 



FIG. 2: (Color online) Evolution of the population on the 
k — quasi-momentum state {F2{k = 0)) after the sec- 
ond expansion. The parameters are {N , Mq , Mi , M2) = 
(50,100,200,300). The dotted line indicates the result with 
the "initial" state being the (first) generalized Gibbs ensem- 
ble (GGE) density matrix Hggei. Other lines correspond to 
results with the "initial" states being tp{ti), with the value 
of ti varied. The markers on the right ends of the lines in- 
dicate the predicted values of the second GGEs. The insert 
shows the time evolution of the population on the k = 
quasi-momentum state {Fi{k — 0)) after the first expansion. 



(to the second expansion) tpiti) replaced by ^ggei- Note 
that since the latter is time independent, this necessarily 
requires that the subsequent dynamics be insensitive to 
the specific value of ti as long as it is large enough to 
belong to the steady regime. 

In the intervals oi t < 0, < t < ti, and t > ti, the 
volume (number of sites) of the system is Mg, Mi, and 
M2 , and thus the corresponding Hamiltonians will be de- 
noted as Hq, Hi, and H2, respectively. They are of the 
form H, = -JEfrLlib]bj+i + 6]+ifej), < z < 2. Here 
J is the hopping strength, and Li = {M2 — Mi)/2 + 1 
and Ri = {M2 + Mi)/2 denote the left- and right-most 
sites accessible to the bosons, respectively. The creation 
and annihilation operators satisfy the usual bosonic com- 
mutation relations plus the hard-core constraint 6| = 

b^^ = 0, so that each site can be occupied by at most 
one boson. By using the Jordan- Wigner transformation 

{b],bj) = {c],Cj)Y{^jTlie~''"''i''''' , where c-, (c]) is the 
fermionic annihilation (creation) operator. Hi is mapped 
to a free fermion one. Hi = —J^j=Li{c\cj+i + c^j^iCj). 
This Hamiltonian can be readily diagonalized as Hi — 
efr{{^^^, with = -2JC0S {^q/{M, + 1)) and 
^ = VSTEf=L.c,sin(g7r(j - + 1)/(M, + 1)). 
The initial state is then simply a Fermi-sea state i/'o ~ 

nr^,.rio). 

From Hq to Hi, the wave function evolves as V'(ii) = 
^-iHiti/h^^^ and the (first) GGE density matrix is de- 
fined as Eggei = ■^exp[-X;*=i A^^^4^^'''4^^]. Here 
the parameter Aq^'* is determined by the initial state. 
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tr{c, 



(l)tjl)77 



(V'olcg^^^Cg^^ l^/io), and 61 is a normal- 



ization factor (or partition function). It is found in [J], 
argued in 12|, and verified in Fig.[2]below that for ti large 
enough, the momentum distribution (or populations on 
the quasi-momentum states, here k — — Mi/2, — Mi/2 + 
I,-- - ,Afi/2-l) 



Fi{k) 



1 



Ri 

E 



^-i27rk{j-j')/Mi,-\ 



(9) 



with respect to il>{ti) can be accurately reproduced by 
using Sggei, i.e., (i/'(ti)| A (fc)|i/'(ii)) ^ tr(Fi (fc)Sggei)- 

Now from i/i to H^, the 772-evolved wave function at 
t = ti+t2 isgivenby 7/)(ii+i2) = e-'"^*^^'^'il;{ti). For our 
purpose, we replace the "initial" state ip{ti) by '^ggei and 
define the i72-evolved GGE density matrix E!ggei(i2) = 
e"'^=*2/^Sggeie'-^2*2/'*. We then study the momentum 
distribution (fc = -M2/2, -A/2/2 + 1, • • • , M2/2 - 1) 



, R2 



-i27rfcO-i')/A^2^t^^^. 



(10) 



with respect to ip{ti + 12) and Sggei(i2). The results are 
shown in Fig. [2j 

In the insert of Fig. [21 we see that after the first expan- 
sion, the population on the fc = quasi-momentum state 
{tp{ti)\Fi{k = O)l^(ii)) relaxes to the steady value pre- 
dicted by the GGE density matrix Sg^ei eventually. This 
proves the predictive power of the GGE after the first 
expansion. What Fig. [2] highlights is that, if the time of 
the second expansion ti is chosen to belong to the steady 
regime, the later evolution of the population on the k ~ Q 
quasi-momentum state {i^{ti + t2)\F2{k = Q)\il;{ti + 12)) 
can be accurately reproduced by tr{F2{k = 0)Eggei{t2))- 
Their lines coincide with each other not only in the 
asymptotic limit but even on details during the tran- 
sitory period. Note that since the latter is independent 
of ti, this necessarily implies that the former is insensi- 
tive to the value of ti, as is indeed the case. Overall, 
Fig. [2] is a remarkable demonstration of the fact that the 
GGE density matrix Sggei shares with the relaxed state 
ipiti) not only the value of the momentum distribution, 
but also the response to a second quench. Or in the per- 
spective of the state ?/'(ii), it has relaxed to be virtually 
indistinguishable from the GGE state Sggei, neither by 
static nor dynamical criterions. 

In Fig. [21 we have also studied whether the steady 
value of F2{k — 0) after the second quench can be 
described by a second GGE density matrix Egge2, 
which is defined as Sgge2 = ^ exp[- J2q=i Ag^^c^^^^Cg^^], 



(2) 

with the parameter Xq determined by the con- 
dition tr(4^-'^4^^Sgge2) = if (4^^^4^'Sggel) Or 

tr(cf ^4'^2gge2) = (i/^(ti)|4'^^4'V(ii)) depending 
on whether the "initial" state is ^ggei or ilj{ti). The 
result is that the second GGEs do predict the steady 



values correctly; moreover, they agree with each other 
very well. This is one more evidence that the relaxed 
wave function '(/'(ii) is virtually indistinguishable from 

the GGE Sggel. 

In summary, we have investigated and verified the rel- 
evance of the GGEs in the dynamical response of the 
two integrable models of transverse Ising model and one- 
dimensional hard-core bosons. Once having relaxed to 
have its properties correctly predicted by the GGE, the 
system behaves as if it were indeed in the GGE state — 
its response to the second quench can be accurately re- 
produced by the GGE even to details. On one hand, 
this result is a welcome complement to previously es- 
tablished result that the GGEs are relevant in predicting 
the static properties of the systems after the first quench. 
The two now combine to present a more complete story 
of the GGE and beckon more confidence on it. On the 
other hand, this result also gives us a sense of "dynam- 
ical typicality" [31, which is also observed in the (non- 
integrable) Bose-Hubbard model previously IJ] . Finally, 
though here we have been dealing with integrable sys- 
tems only, a lesson may also be drawn for non-integrable 
systems. A closed non-integrable system might well be a 
pure state yet virtually indistinguishable neither by static 
nor by dynamic criterions from a canonical ensemble. 
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